1. IntroductionGraphene plasmon is attracting increased research attention because of their potential applications in technology, especially in optoelectronics.[1–3] Many theoretical studies have focused on terahertz graphene plasmon[4–6] because it exhibits many unusual behaviors, such as high wave velocity and long propagation distance.[4,7,8] Engineering a tunable graphene plasmon device is one of the key experiment tasks to determine their superior characteristics. Tunable terahertz plasmon in graphene ribbons/resonant electromagnetic (EM) cavities has recently been investigated to develop plasmon-based terahertz devices. However, the coupling between graphene and incident terahertz wave is relatively weak.[9,10]
The adoption of metallic grating and resonant cavity is one of the best solutions, and such a scheme has been studied in our previous paper.[6,11] The upper grating serves as a superior coupler to excite graphene plasmon.[6,11,12] The interaction between graphene plasmon and EM wave can be enhanced by a Fabry–Pérot cavity.[6,11,13] The relation between graphene plasmon frequency and the grating width has also been quantitatively analyzed.[11]
This study investigates the voltage-tunable terahertz plasmon in a grating-coupled two-dimensional electron density (2DEG) embedded in a Fabry–Pérot cavity, using the terahertz time-domain spectroscopy (THz-TDS). Transmission theoretical calculation on the graphene plasmon agrees well with the experimental results. Methods to enhance the tunability of the prepared device are discussed in the next sections.
2. ExperimentsThree main building blocks were involved in this device, a graphene sheet, a Fabry–Pérot cavity, and a grating, as shown in Fig. 1(a). The graphene was synthesized by chemical vapor deposition (CVD) method under low pressure and high temperature. After that, it was transferred onto a silicon substrate (THz permittivity, ε1 = 11.9) with a 300-nm oxide layer (THz permittivity, ε2 = 4) by using the PMMA transfer method. Raman measurements were performed in the as-prepared sample. Results show that the transferred graphene is single layered with high quality[14,15] (see Fig. 1(b)). A thin film of Al2O3 was deposited on the sample by atomic layer deposition (ALD)[16] subsequently, whose chamber temperature was set to 140 °C and 280 rounds of deposition were carried out. The thickness of Al2O3 film deposited on graphene was about 25 nm. These three elements (Al2O3/SiO2/Si) introduced above make up the main body of the cavity. Next, the 100-nm thick gold grating gate, coupling the terahertz EM wave to the grapheme plasmon and regulating the electronic density of graphene, was fabricated on top of the Al2O3 layer. The grating constant L was 8.7 μm when the length of the individual gold finger W was 4.7 μm (see Fig. 1(c)). The source and drain electrodes were produced simultaneously when the grating gate was prepared. Besides, the device was encapsulated using a chip carrier with 24 pins. The overall thickness and the total functional area of the sample were D1 = 190 μm and S1 = 4 × 4 mm2, respectively. Finally, the THz-TDS technique allowed us to obtain the transmission spectrum at different gate voltages. The incident light polarized along the grid fingers, as shown in Fig. 1(a).
3. Results and discussionTransmission of the bare cavity and the grating-coupled graphene plasmon device are shown in Fig. 2(a). The bare cavity is an insulator with constant transmission level in the THz region. However, the transmission spectra measured on the graphene device with resonant cavity shows that a prominent peak locates at 1.511 THz. This phenomenon implies that THz absorption is observably affected by graphene. Thus, extinction spectra (1 − T/T0) for grating devices with VTG = −0.8 V, 0 V, 0.8 V were then investigated.
To compare with the experimental data, we calculate (1 − T/T0) of the graphene device theoretically by[11]
where
T0 and
T are the transmissions passing through the bare cavity and the sample with the graphene, respectively.
Z0 is the vacuum impedance, and
ns is the effective refractive index of the substrate. The solid curves in Fig.
2(b) are fitting results according to Eq. (
1) and the Lorentz conductivity. We use Lorentz conductivity for the graphene layer as
[10,11]
where
ω is the radian frequency of the incident terahertz light,
n is the carrier density, and
Γp is the plasmon resonance width. The resonance frequency is given by
[10,11]
where
ε and
kG are effective dielectric constant and plasmon wave vector, respectively. The
σuni = π
e2/2
h is the well-known optical conductivity, and the Femi energy is given by
EF =
hvF(
πn)
1/2,
vF = 1 × 10
6 m/s. If the sample is irradiated vertically by the light, the plasmon wave vector is determined by
kG = m
π/
W with
m = 1,2,3,…
[11,17]The corresponding solid curves in Fig. 2(b) are fitting results according to Eqs. (1)–(3) with n and Γp as two fitting parameters. Three prominent features for the data are presented in Fig. 2(b). First, the plasmon resonance width of the as-prepared sample is approximately 3.5 THz, whereas the carrier density n ranges from 7.16 × 1012 cm−2 to 8.36 × 1012 cm−2. Second, an obvious blue shift of the plasmon resonance frequency with decreasing top-gate voltage is observed: resonance frequencies are located at 1.523 THz, 1.511 THz, and 1.506 THz for voltage of 0.8 V, 0 V, and −0.8 V, respectively. Using a larger range of effective top-gate voltage, terahertz resonance in 1–2 THz can be straightforwardly achieved (see Fig. 5). Third, the light–plasmon coupling in graphene is remarkably strong; transmission extinction at the plasmon resonance for VTG = 0 V is more than 80%. A higher extinction can be achieved readily by enlarging the top-gate voltage VTG.
The electrical properties of the engineered graphene plasmon device are discussed by performing a fullest direct current (DC) measurement. Figure 3 shows the equivalent current model in the DC measurement. First, the output character is measured to obtain the device resistance (see Figs. 3(a) and 4(a)). The channel current of the graphene device is found to be linear with the output voltage VDS, and the output resistance is proved to be a constant of 1639 Ω. Second, the leakage character of the fabricated graphene plasmon device is obtained in Figs. 3(b) and 4(b). In the measurement, the source and drain electrodes are linked to the ground. Figure 3(b) shows the leakage measurement equivalent model. The top-gate structure can be regarded as a parallel between the top-gate capacitance Cx and junction resistance RG because of the top-gate leakage. In DC measurement, the current can only flow through the resistance RG. The leakage current results IG1 and IG2 show an almost perfect symmetry to the origin of coordinates (Fig. 4(b)). The total leakage current is IG = IG1 + IG2 ≈ 1.5 μA when VTG = 0.85 V. Thus, the junction resistance RG is approximately 0.6 MΩ, far less than that of small-size top-gate media. Meanwhile, the value of IG1 is almost the same as that of IG2. Assuming two graphene channel resistances close to the source and drain electrodes are equal, RS ≈ RD.
Third, our double DC measurement of the prepared graphene device for different top-gate voltages will explicitly indicate the transfer character, as shown in Figs. 3(c), 3(d), and 4(c). From the values of ID, IS, and IG extracted from the DC measurement (inset in Fig. 4(c)), the actual top-gate voltage VG set on the capacitance Cx at different gate voltages can be obtained in Fig. 4(c). Here, VG = VTG − V (see Figs. 3(c) and 3(d)), with graphene channel voltage V = IDRD + RC/2) = VDS − IS(RS + RC/2). To further simplify, RS + RC/2 = RD + RC/2 = VDS/(IS + ID), V = VDSID/(IS + ID). As a result, VG = VTG − VDSID/(IS + ID), and the actual channel current due to the graphene carrier can be expressed by IDS = V0/Rt, where the whole device resistance Rt can be given by Rt = RD + RS + RC = 2VDS/(IS + ID), and V0 is set as 10 mV. The theoretical fitting for our transfer characteristic data will be discussed below. The carrier concentration (electrons or holes) in the graphene channel regions n related to the residual carrier concentration n0 and the Dirac-point voltage VDirac can be approximated by , ,[18] where Cx and n0 can be set as 122.4 nF/cm2 and 1 × 1012 cm−2 in our fitting model, respectively. The total device resistance Rt is given by Rt = Rcontact + L /(Weμn),[18] where the contact resistance is Rcontact = RS + RD and μ represents the mobility. Finally, the channel current is given by IDS = V0/Rt. The measured IDS against VG (symbols) along with the fitting results (solid lines) is shown in Fig. 4(c). Then, the relevant parameters can be extracted: VDirac = −10.5 V, μ = 1500 cm2⋅V−1⋅s−1 and Rcontact = 770 Ω.
However, the variation range of the graphene plasmon frequency is very small, as shown in Fig. 5. The calculating parameters were set as the same as those in Fig. 4(c). In our experiments, tiny top-gate voltage range was applied to the device. Otherwise, when the larger top-gate voltage is applied to modulate the carrier density of the graphene channel, the leakage current can reach up more than 1.5 μA (see Fig. 4(b)). Therefore, the top-gate voltage can only be set as small as possible to avoid the large leakage current to breakdown the top-gate capacitance. Two aspects should be considered to develop a device with a larger plasmon frequency range. First, defects of the top-gate media can be decreased by enlarging the growth temperature of ALD and applying the graphene prepared using molecular beam epitaxy (MBE).[19] As a result, the sustaining capability of the top-gate media can be improved directly. Second, the thickness and permittivity of top-gate dielectric layers can be optimized to enlarge its capacitance. Consequently, the range of graphene plasmon frequency can also be expanded (see the inset of Fig. 5).